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equation depending on the effective mass and a suitably chosen relativistic potential, one obtains a 
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I. INTRODUCTION 

The recent advances in nanostructure technology al- 
low the fabrication of electronic devices which inherently 
experience quantum effects in their operation. There is 
a need for theoretical tools, not only to understand the 
behavior of actual nanoscale devices, but also to initiate 
future research and developments. 

The simplest non-relativistic quantum modeling of 
nanoscale semiconductor devices is based on the 
Schrodinger equation written in solid state domains 
where the potential is constant and the influence of the 
lattice is encapsulated in the value of effective electron 
mass. When the devices are made from semiconductor 
heterostructures, there are many such domains separated 
among themselves by interfaces where besides the step in 
potential we also have to consider the discontinuity in ef- 
fective electron mass In the '60 the problem has been 
solved by applying appropriate boundary conditions ^ , 
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where ip^, and ^^j, m^j are the electron wave func- 
tion and effective mass of the electron to the left {L) 
and right (i?) side of a given interface. The results ob- 
tained with the above-described procedure are in agree- 
ment with experimental data and the method is widely 
used in electron transport computations through semi- 
conductor heterostructures The advantage of this 
model is that the transmission coefficient can be calcu- 
lated using the simple and elegant method of the transfer 
matrix 0, However, the assumed boundary condi- 
tions are imposed somewhat artificially in order to 
conserve the particle current without to have a deeper 
physical motivation. 

Another attitude is to start with the Dirac equation 
even though it is clear that the relativistic effects have 
to be very small. Nevertheless, the relativistic linear 
dependence between energy and mass could offer some 
technical advantages for finding appropriate connection 
conditions at interfaces where the potential and the ef- 
fective mass present discontinuities. Few years ago, the 
one-dimensional Dirac equation was successfully used for 



treating problems with variable mass avoiding several dif- 
ficulties of the non-relativistic theory 6]. 

In this paper we would like to continue this study us- 
ing the normalized plane wave solutions of the three- 
dimensional Dirac equation in helicity basis. Our pur- 
pose is to derive the relativistic version of the transfer 
matrix method for the motion in a fixed direction of a 
Dirac electron with point-dependent effective mass, pass- 
ing through rectangular barriers of arbitrary profile. We 
show that the use of the Dirac equation allows one to im- 
pose simple connection prescriptions at interfaces. How- 
ever, the price to pay for working with variable mass is 
that there are many energy scales corresponding to dif- 
ferent mass values. For this reason we need to rescale the 
experimental potential if we want to measure the ener- 
gies with respect to an unique energy scale. The rescaled 
potential will be considered the appropriate relativistic 
potential of our problems. We show that only in this 
way the non-relativistic limit of our approach recovers 
the results derived from Schrodinger equation with the 
conditions 1^. 

The paper is organized as follows. In the second section 
we present the well-known plane wave solutions of the 
Dirac equation in the helicity basis. The next section is 
devoted to the problem of one-dimensional rectangular 
barriers of any shape allowing us to find the relativistic 
transfer matrix in section four. Finally, it is shown that 
in the non-relativistic limit this matrix becomes just the 
desired traditional one lal. 



II. PLANE WAVES 

Let us consider the Minkowski space-time in a frame 
of coordinates i^, ... — 0, 1, 2, 3) and the metric -q = 

diag(l, — 1, — 1, — 1). In natural units (with fi = c — 1) 
the time is x^ = t while the space coordinates, x^ = x, 
x^ = y and x^ = z, are the components of the vector 
X. In this frame, the relativistic quantum motion of an 
electron of mass m and charge — e, in an arbitrary ex- 
ternal electromagnetic field A^, is governed by the Dirac 
equation 0, 

7^(za^-eA^)V'-mV' = 0, (2) 
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that produces the conserved current (in units of — e) 

f = ^j^^P , (3) 

where ip = ifj^^'^ is the Dirac adjoint of the spinor ip. 
In what foUows, we take the 7- matrices in the standard 
representation (with diagonal 7°) 0- 

Here we are interested to study the quantum modes in 
the particular case of a space domain D where A{x) = 
and eAo{x) = V = const, for any x G D. In this do- 
main the Dirac equation can be analytically solved and 
different quantum modes can be well-defined using com- 
plete sets of commuting operators. Thus the plane wave 
solutions are eigenspinors of the complete set of commut- 
ing operators {Ed,P,W} constituted by the Dirac op- 
erator, Ed = ij^df_i — j'^V, momentum P = iV, and 
the Pauli-Lubanski operator W = 2P ■ S. The cor- 
responding eigenvalues, m, k and A, define the plane 
wave spinor of positive frequency, momentum fc, energy 
E{k) = Vm^+k^ +V and hehcity A that reads ^Sl 
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E{k) -V + m^x{k) 



X e 



E{k) - V - m ^x{k) 

— iE{k)t+ik-x ^^■j 



We denoted by £,\ {k) the normalized Pauli spinors of the 
helicity basis that satisfy k ■ a £,x{k) — A|A:|^A(fc) and 
[£,\{k)]~^^\'{k) — (5a, A' (where Ui are the Pauli matrices 
and A = ±1). One can verify that each solution is 
normalized as tp^ ^tjjj: y — Sx^y and produces the current 
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FIG. 1: A sequence of potential steps. 



while in the domains Din and Dout its mass is just the 
bare mass m. 

In special relativity the energy scale depends on the 
value of the rest mass while the electromagnetic poten- 
tial is defined up to a gauge. Therefore, in problems 
where this mass is replaced by a point-dependent effec- 
tive mass, we could introduce an unique energy scale only 
by choosing suitable gauge fixings, dealing with the dif- 
ferent values of the effective mass. In these conditions we 
are encouraged to consider in each domain Di the rela- 
tivistic potential Vi instead of the experimental one Vi. 
The relation among these potentials has to be derived 
from a natural supplemental condition which will fixe up 
the gauge in the domains Di. 

In any domain Di there exists a plane wave solution of 
energy E and helicity A propagating in the sense of the 
positive semiaxis z. 



1 



^ — iEt~\-ikiZ 



(6) 



1 — ^ Ifcl 

^■ = ^^^,A(fc-7)^J,A = - 



(5) 



along the direction k. 



III. ONE-DIMENSIONAL MOTION 

The general results presented above help us to write 
down the solutions of simpler one-dimensional problems 
along the third axis. Of a special interest is the problem 
of the electron moving through a system of N rectangular 
barriers of arbitrary shape. In general, the system of bar- 
riers is constituted by N domains Di = [zi,Zi_|_i] where 
the potential V{z) takes constant values Vi (Fi^). These 
domains are limited by plane interfaces at fixed points, 
Zi, Z2, ^Af+i, among them those from zi and zat+i rep- 
resent the interfaces between the system of barriers and 
the domains outside, denoted by Din = Dq — (— 00, zi] 
and, respectively, Dout = Dn+i — [-2^+1,00). It is natu- 
ral to consider that in these latter domains the potential 
vanishes, Vin = Vout = 0. In addition, we assume that 
in each domain Di the electron has the effective mass rrii 



which depends on the constants fcl^-* = y/pf- 
and scalar momentum 



— — 



{E^Vi)^-mj 
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(7) 



We note that in this case the helicity spinors coincide 
to those of the spin basis since the spin is projected on 
the third axis. Consequently, the two-component spinors 
^A take the usual form ^1 = (1,0)-^ and = (0,1)-^. 
The plane wave solution with the same E and A but 
propagating in the opposite sense reads 
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The conclusion is that, in a domain 1?^, the most general 
plane wave solutions of energy E and helicity A are given 
by the linear combinations 

where Ai and Bi are arbitrary complex numbers. Each 
solution © gives the total current 



(10) 
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which does not depend on heHcity. 

Finally we can establish the relation among the rela- 
tivistic and experimental potentials assuming that in a 
domain Di the momentum ki vanishes only when the to- 
tal non-rclativistic energy Enr — E — m, calculated with 
respect to the bare mass m, equals the experimental po- 
tential Vi. Therefore, according to Eq. O we obtain the 
form of our relativistic potentials 



Vi = Vi + Srui , 



(11) 



where Sruj 



The global solution of energy E and helicity A is contin- 
uous in each point Zi which means that we must impose 
the conditions 



(17) 



for i=l,2,...,Af-|-l. After a few manipulation we find 
that these conditions lead to simple relations among the 
associated vectors, 



where the matrices 



i = 1,2, 



,N+1. 



(18) 



IV. THE TRANSFER MATRIX 

In what follows we shall derive the transfer matrix in 
the pure scattering case without wells or tunneling ef- 
fects. This means that the energy satisfies the condition 
E > Eq ^ sn^{Vi m I z = 0, 1, N + I) and h take 
only real values. In addition, we specify that the global 
solutions we consider here have the same fixed energy E 
and helicity A in all the domains Di. This correspond to 
the experimental situation when the interactions able to 
produce spin-flip are absent. 

In problems involving many domains Di it is difficult 
to manipulate solutions of the form Q. For this reason 
we replace_these solution by associated two-dimensional 
vectors 



Vi{z) 



Bie-'^^ 



(12) 



which carry all the information we need for calculating 
the currents Hll)|) . Indeed, we observe that these currents 
can be expressed only in terms of Vi{z) as 
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[vi{z)\^ a^Viiz) 
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(13) 



Thus the vectors H12|l become the basic elements of the 
relativistic formalism of the transfer matrix for rectan- 
gular barriers 5.]. In the domains Dq and Djq+i, where 
the potential vanishes and the mass is to, the spinors 
and "^^^ have the general form given by Eqs. ((HI, 
and with the same momentum, fco = fcAr+i = 

'^^^ - VE ±m. 



m?, and constants, fcp^' = fcj^^]^ 



k = 

These spinors are associated to the vectors 



Vin{z) = Vq{z) = 



A ■ p^- 

Din^ 



-ikz 



Vout[z) = VN+l[z) = ik^ 



(14) 
(15) 



Now, the problem is to find the transfer matrix, M, which 
transforms the out vector into the in one as 

V(i{zi) = MvN+i{zN+i) , (16) 
allowing one to calculate the transmission coefficient. 
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depend on the constants 
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The last step is to introduce the translation matrices 



T — 



-iki{zi+i-Zi) Q 

giki{zi + i-Zi) 



(21) 



which transform Vi{zi) into Vi{zi) = TiVi{zi+i). With 
these elements we can write down the final expression of 
the relativistic transition matrix 
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M 
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(22) 



Let us observe that for E > Eq, when ki are real num- 
bers, the matrices Mi = defined by Eq. H18|l have 
the property 



TOt_i ki 
ki-i TOj 



(73 



(23) 



which guarantees the conservation of the total current, 
jm = ji = ■•• = ji = •■■ = jout, calculated according to 
Eq. (|13|) . In these circumstances, taking Bout = we 
have l^inp — |i?mP = l^outP which allows us to define 
the transmission coefficient 



T = 



\Mn 



(24) 



We note that T results to be a function only of energy, 
being independent on the helicity of the electron passing 
through the rectangular barriers. 

Of course, the function T{E) calculated here is defined 
only on the domain E > Eq. However, starting with the 
present theory, the extension to energies smaller than Eq 
can be done but this requires specific treatment because 
of the wells producing discrete energy levels or tunneling 
effects which need to be treated with specific methods. 
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V. CONCLUSIONS 

Here we constructed the relativistic version of the 
transfer matrix for the Dirac electron moving through 
rectangular barriers, in a similar manner as in the non- 
relativistic theory based on the Schrodinger equation . 
Our approach allows one to calculate the transfer matri- 
ces using the same rules but with matrices Mi of different 
forms. 

Let us see what happens with our theory in the non- 
relativistic limit, for small values Enr — Vi<^mi. In this 
limit the quantities fc,- ^ = V^nr ~ Vi remain unchanged 
but we have fc,-^'' — > y/2mi. Consequently, we find that 
Tj-"*^' — > 1 remaining with the terms 

which coincide to those of Refs. (al at least in the domain 



E > Eq considered here. Thus, the general conclusion is 
that the non-relativistic limit of our approach based on 
the three-dimensional Dirac equation with the relativistic 
potentials reproduces identically the results of the 
traditional method based on the Schrodinger equation 
and conditions Q|. Moreover, relativistic corrections can 
be also calculated but so far these seem to be small in 
the usual regime the electronic devices work. 

In other respects, the results obtained here indicate 
that the use of the Dirac equation could be helpful in 
other problems concerning the motion of electrons in 
semiconductor heterostructures as suggested in Ref. 0. 
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